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1 — Abstract 

t-H ■ 

Natural cilia are hair-like microtubule-based structures that are able to move fluid at low Reynolds 

00 '. 

number through asymmetric motion. In this paper we follow a biomimetic approach to design 

m ; 

artificial cilia lining the inner surface of microfluidic channels with the goal to propel fluid. The 

o : 

ON . artificial cilia consist of polymer films filled with magnetic nanoparticles. The asymmetric, non- 
O ■ 

reciprocating motion is generated by tuning an external magnetic field. To obtain the magnetic 

field and associated magnetization local to the cilia we solve the Maxwell equations, from which the 

magnetic torques can be deduced. To obtain the ciliary motion we solve the dynamic equations of 

motion which are then fully coupled to the fluid dynamic equations that describe fluid flow around 

the cilia. By doing so we show that by properly tuning the applied magnetic field, asymmetric 

ciliary motion can be generated that is able to propel fluid in a microchannel. The results are 

presented in terms of three dimensionless parameters that fully delineate the asymmetry and cycle 

time as a function of the relative contribution of elastic, inertial, magnetic and viscous fluid forces. 
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1. Introduction 

In the biomedical field there is an increasing need for miniaturized lab-on-a-chip analysis sys- 
tems that are able to analyze small quantities of biological samples such as biofluids (e.g. blood, 
saliva, urine) 
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ly, yj. These so-called biosensors are microfabricated total-analysis-systems that 
typically consist of a system of microscopic channels, connecting microchambers (the labs) where 
dedicated tests are carried out. Classical means for fluid-propulsion do no longer suffice at these 
small length scales, which has led to a search for new methods dedicated for fluid propulsion at 



the micron-scale, such as micropumps 



actuation, as in electro-osmotic 



[lOj, syringe pumps 

mm 

or by exploiting electro-magnetic 



la. 14] and magnetohydrodynamic devices [l|, [6|. The use of elec 



Diting 

0l 



trie fields, however, in transporting biological fluids (which usually have high conductivity) may 



In 



induce heating, bubble formation and pH gradients from electrochemical reactions 
this work, we explore a new way to manipulate fluids in microfluidic systems, inspired by nature, 
through the magnetic actuation of artifical cilia. 

Akin to size-effects in the mechanical properties of solid materials, the dependence of size in 
fluid dynamics is captured by the Reynolds number. At small length scales and low Reynolds 
number, gravity does not play a role and fluid dynamics is usually dominated by viscosity rather 
than inertia. An important consequence of this is that fluids can only be propelled by motions that 
are initiated by actuators whose movement is cyclic but asymmetric in time. Nature has solved this 
problem by means of hair-like structures, called cilia, whose beating pattern is non-reciprocating 
and consists of an effective and a recovery stroke, as shown in Fig. [TJ In this work we design artificial 
cilia that can be actuated by an external magnetic field. The artificial cilia are thin films consisting 
of a polymer matrix filled with magnetic nanoparticles. Depending on the nature of the particles, 
the film can be super-paramagnetic or ferromagnetic with a remanent magnetization. The applied 
magnetic field is uniform but its magnitude and direction can be manipulated in time to get the 
desired asymmetric motion. The Lagrangian model for the cilia is based on a dynamic finite- 
element representation, accounting for elastic, inertia and drag forces in a non-linear geometry 
setting. Simultaneously, Maxwells equations are solved at each configuration, yielding the local 
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magnetic field and magnetic induction from which the magnetic forces can be obtained. This 
magneto-mechanical model allows studying the overall deformation of the cilia actuated by an 
external magnetic field as a function of the mechanical and magnetic material parameters of the 
films microstructure. The Lagrangian solid model is also coupled to an Eulerian formulation of the 
fluid, which enables to generate fluid flow through the magnetically-induced film motion. The goal 
of this paper is understand the multiphysics interplay between magnetostatics, solid mechanics and 
fluid dynamics and to explore ways to exploit this interplay to design non-reciprocating motion in 
artifical cilia. 

The paper is organized as follows. In section 2 we discuss the magnetomechanical model 
consisting of a finite element discretization of the dynamic principle of virtual work that is explicitly 
coupled to a discretized form of an integral formulation of Maxwells equations. In section 3.1 and 
3.2 we study several fundamental magnetic loading situations in which the interaction of elastic, 
inertial, magnetic and viscous drag forces is summarized in a deformation mechanism map that 
features three distinct regimes of operation. In section 3.3 we identify four different configurations 
that show asymmetric motion, followed by a parametric study in section 4 to explore the efficiency 
of ciliary motion as a function of the three governing dimensionless numbers. Finally, we couple the 
magnetostatic model a fluid-dynamics model and show that fluid can be propelled and that a linear 
relation exists between the asymmetry of the ciliary motion and the amount of fluid propelled per 
cycle. Section 5 contains the conclusions. 

2. Method 

The proposed actuators are polymers films filled with magnetic particles. Depending on the 
magnetic nature of the particles, the film can either be super-paramagnetic (SPM) or permanently 
magnetic. The thickness of the film is much smaller than its length, which allows us to use Euler- 
Bernoulli kinematics. In the following two Sections we will discuss the numerical methods used to 
solve the equations of motion for the dynamic deformation of the film and Maxwell's equations for 
the magnetostatics problem. 
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Effective Stroke 




Figure 1: Asymmetric motion of a cilium. Instances 1-4 show the effective stroke and 5-8 show the recovery stroke. 
The asymmetry in the motion of cilia during forward and return stroke propels the fluid effectively in one direction. 

2.1. Equations of motion 

As a starting point for the Euler-Bernoulli beam element formulation we use the principle of 



virtual work 



In this weak form of the equations of motion, the virtual work of the external 



forces is equal to the internal work and is given by 



8Wt 



5Wt 



(1) 



with 



swl 



hit 



(a5e + p{u5u + vSv)) dV, 



(2) 



where u and v are the axial and transverse displacements along the beam length (axial coordinate 
x) and p is the density of the film. Furthermore, a is the axial stress and e is the corresponding 
strain, given by 



du 1 f dv 
dx 2 \ dx 



d 2 v 
dx 2 



Its first variation is thus 



. d5u dv dSv d 2 5v 
oe = — h - — y -z—z = oe 



dx dx dx 



dx 2 



yb~x- 



(3) 
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By substituting the strains and defining J odA = P and — J oydA = M (A is the area of the cross 
section), the internal virtual work at time t can be written as the sum of an elastic and an inertial 
part 



SW^ t = J [P^t 1 + M^x 1 + pAivtdu 1 + v t 5v t )Jdx. (4) 
The internal virtual work at time t + At is written as 

SW*+ At = J ( P t+^Se t+At + M t+At 5 X t+At + pA(u t+At 5u t+At + v t+At 5v t+At )) dx. (5) 
The corresponding external virtual work is 



8W* M = j ^f x +At 5u t+At + f y +At 6v t+At + N^^^j Adx 
+ J (ti +At 5u t+At + tl +At 5v t+At ) bdx, 



(6) 



where f x and f y are the magnetic body forces in axial and transverse directions, N z is the magnetic 
body couple in the out-of-plane direction, t x and t y are the surface tractions due to the fluid drag 
and b is the out-of-plane thickness of the film. We now expand the elastic part of the internal work 
linearly in time by substituting (Q t+At = Q 1 + AQ) for the field parameters in Eqn. [5] which gives 

5Wl+ At = J [(P'Se 1 + M'Sx 1 ) + (AP5e* + AMSx 1 ) + P'ASe + pA(m* +a '5m* + i/ +A W)] d(ff,) 

in which terms of order higher than one are neglected. The axial and transverse displacements 
are linearly and cubically interpolated in terms of the nodal degrees of freedom, the displacements 
and rotations, 

u = N u p, v = N v p, (8) 

where 

p = {«i Vi <jMo U 2 V 2 4>2lo} T , 
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N u and N v being the standard interpolation matrices, given in Appendix El By using the following 
notation 

du _ dv _ <9 2 m _ 

— = — = — = C>, 9 

ax ax ox 2 

and the constitutive relations 

AP = EAAe, AM = EIA X , (10) 



with E being the elastic modulus and I being the second moment of area defined as I — j^bh 3 , 
the internal virtual work can be written as 



Sp T J [P* (Bj + B v T B v p) + M l Cj] dx 
+ 5p T J [EA (Bj + BjB vP ) (B u + B v B v p)] dxAp 
+ 5p T J EICjC v dxAp + 5p T J P l B v T B v dxAp 
+ 5p T J P A{NjN u + NjNJdxp* 



(11) 



■J+At 



By choosing the domain of integration to be the current configuration (i.e. using an updated 
Lagrangian framework), the total displacements are zero, p — 0, and we get 

5Wf+ At = 5p T fl t + 5p T KAp + 5p T Mp t+A \ (12) 

where 

fl t = J [P l Bj + M t Cj]dx 

is the nodal internal force vector, 



K = I EAB u r B u dx + / EIC v r C v dx + / P l B v T B v dx 
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is the stiffness matrix, the first two terms of which represent the material stiffness and the third 
term represents the geometric stiffness, and 



M = J P A(NjN u + N v T N v )dx 



is the mass matrix. The external virtual work is 



5W£ At =5p T J [(f x +At Nj + f y +At N v T + N t z +At B v T )A + b (t^Nj + ^Nj)] dx, 



(13) 



-s P T fl At - 



By equating the internal and external virtual work and noting that the resulting equation holds 
for arbitrary Sp we get 



KAp + Mp^ = f c t Ai -f mt . 



(14) 



The motion of the film with time is obtained by solving Eqn. [TH with appropriate initial and 
boundary conditions. Newmark's method is employed to integrate Eqn. [TH in time, viz, 



P 



t+At 



p" + p*Ai + -At 2 [(1 - 20)p x + 2pp t+At ] 



p t+At = p* + At[(l-7)p t + 7 p* + - At ] 



(15a) 
(15b) 



where 7 and (3 are integration parameters. By using Eqn. I15a[ Eqn. [TH can be written as 



(At 2 (3K + M)p t+At = f c t At - /L " 



p* + -At(l - 2(3)p t 



(16) 



the solution of which will give the nodal accelerations, which are integrated using Eqn. I15al and 
Il5bl to get the nodal displacements and velocities. 

The drag forces of the fluid on the film are accounted for through the surface tractions in 
Eqn. [131 The tractions are proportional to the velocity in the low Reynolds number regime, with 
the normal and tangential coefficients of proportionality (the drag coefficients) denoted by C x and 
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Cy. The tractions are calculated from the velocity at time t. The magnetic forces come in through 
the body forces and couples in Eqn. [131 These are calculated by solving Maxwell's equations on 
the current configuration, i.e. the configuration at time t (see Section 1272]) . 

2.2. Magneto statics 

Maxwell's equations for the magnetostatic problem with no currents are 

V B = (17) 

V x H = 0, (18) 

with the constitutive relation 

B = fi (M + H), (19) 

where B is the magnetic flux density (or magnetic induction), H is the magnetic field, M is the 
magnetization which includes the remnant magnetization, and /xq is the permeability of vacuum. 
Substituting Eqn. [19] into Eqn. [T71 yields 

V H = -V-M. (20) 

As V x H = 0, a scalar potential <ft exists, such that H = — V0. Substituting this in Eqn. [201 yields 
a Poisson equation for <ft, V 2 </> = —V • M. By taking into consideration the effect of discontinuity 
in the medium, the general solution of the Poisson equation can be found [7(, resulting in 

H(x) = --L v / ";- M <*v + .Lv / v ;- M( :'w-. 

Ait J \x — x \ Air J \x — x \ 

where n' is the outward normal to the surface of V. By assuming that the magnetization is uniform 
inside the volume, so V' • M = 0, the volume integral vanishes, and the field is only due to the 
jump of magnetization across the surface, as reflected by the surface integral in Eqn. [21] 

We now discretize the film into a chain of rectangular segments. On the surface of each segment 
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there is a jump of magnetization. M is uniform inside the segment and zero outside the segment. 
The magnetic field in local coordinates (denoted by ") due to the four surfaces of a segment can 
now be calculated at any position (x, y) by evaluating the surface integral in Eqn. [2TJ resulting in 



tan 



+ tan 



2tt 
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M x ( — tan 



-i 



+ tan 



My |ln 



| In 



2n 

2 



~M y ( \ In 




H-£) 2 + (l + ^) 2 ]-^ ln [^- £ ) 2 + ^ + ^ 2 ]) 



2tt 



M„ 



tan 



-i 



e ; - 



+ tan 



-i 



-M„ 



2tt 



+ 



tan 



-i 



h ; - 

2+^ 



+ tan 



-i 



(22) 



Mx |ln 



H+-) 2 + (-l-^) : 



Mn 



2vr 

2 



2tt 



(f + £) 2 +H-y) 2 ]-|ln [(f + x) 2 +(|-y) 2 ]) 



2tt 

where M x , M y are magnetizations in the tangential (or length) and normal (or thickness) directions, 
h is the thickness and I is the length of the segment. Here, x and y are the local coordinates having 
their origin in the center of the segment. On rearranging, the magnetic field due to segment % can 
be written as 



Gi Mi. 



(23) 



where M = [ M x M y } T ', H = [ H x H y ] T and Gi can be obtained from Eqn. [22j The field due 
to segment i with respect to the global coordinates is 



Hi — R t Hi, 



(24) 
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where Ri is 



cos 9i — sin $i 
sin 9i cos 9i 



(25) 



with 9i the orientation of the segment i with respect to the global coordinates. The field at any 
element j because of the magnetization of all the segments throughout the film is 



N 



Hj = / / + RiGijMi 



(26) 



where Gij properly accounts for the relative positioning of segments i and j, H is the externally 
applied magnetic field, far away from the film, and N is the total number of segments. By rotating 
Hj back to the local coordinates we get 



N 



Hj = RjH + Rj RiGijMi 



(27) 



For the situation of a permanently magnetized film with magnetization Mi, i = 1, . . . , N, Eqn. 
gives the magnetic field in all the segments. 

However, in case of a super-paramagnetic film, the magnetization Mj is not known a-priori, 
but depends on the local magnetic field through 



Mi 



±Hj 

xRjH + ±R] 'RiGijMi 



N 



:28) 



i=i 



X 



with 

Xx o 

Xy 

There are iV similar pairs of equations. In total these are 2xN equations for the 2xN unknown 
magnetizations. This set of equations is solved to get the magnetization with respect to the local 
coordinate frame. Once the magnetization is known the field can be found from Eqn. [271 Then 
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the magnetic flux density can be found by using Eqn. [T9j The magnetic couple per unit volume is 
given by AT" = M x B, which can be obtained in local coordinates through N z = M x B y — M y B x . 
Magnetic body forces due to field gradients will be neglected in the present study. 

3. Results 

3.1. Fundamental loading situations 

In this Section we study the behavior of the film under simple loading conditions. We will con- 
sider a horizontal film that is clamped at the left end and free elsewhere. Three cases are considered: 
(i) a permanently magnetic film with a field applied in the transverse (vertical) direction, (ii) a 
permanently magnetic film subjected to a rotating magnetic field and (iii) a super-paramagnetic 
(SPM) film subjected to a rotating magnetic field. 

The following reference parameters are considered: Elastic modulus E — 1 MPa, drag coeffi- 
cients C x = 10 Ns/m 3 , C y = 20 Ns/m 3 , remnant magnetization for permanently magnetic films 
M = 15 kA/m, susceptibility for the SPM film \ x =4.6, % y =0.8, reference time during which the 
field is applied t re f = 1 ms, length of the film L = 100 /im, thickness h = 2 /im and the magnitude 
of the applied field at t = t re f is -Bomax = HoHo max = 12.5 mT. In the simulations the values of 7 
and P are 1.0 and 0.5, respectively. 



Permanently magnetic film in transverse and rotating field. A magnetic field is applied in 
the vertical (y) direction and is increased linearly from zero to -Bomax in £ = £ re f. Then the applied 



magnetic field is maintained at this value. Fig. 2(a) shows the vertical displacement of the free 
end of the film as a function of time. After time t = £ re f the applied field is kept constant, but due 
to the presence of the fluid and inertia in the film, it reaches a steady state at time t — 1.5 £ re f. In 



Fig. 2(b) the entire film is depicted showing the film deformation in time at five instances identified 



in Fig. 2(a). Next to the film, at the right, the arrow denotes the applied magnetic field B at 
five instances. In addition, the magnetic induction B is plotted by means of arrows along the 



film. Fig. 2(c) shows the variation of the horizontal component B x (with respect to the global 
coordinate axes) along the film. Since no external field is applied in the x direction, the primary 



11 



cause of B x is the permanent magnetization in the film, M. For a permanently magnetic film, 
we can write the magnetic field as the sum of the applied field H and the field generated by the 
magnetization in the rest of the film (the 'self-field'), H = H + H se n, see Eqn. [261 It can be 
clearly seen that at instant 1 when the film is still horizontal, that it is this self-field that tends 
to decrease the magnetic induction near the ends of the film, effectively de-magnetizing the film. 
When the film deforms, the horizontal component of the field, B x , decreases near the free end 



while By increases (see Fig. 2(d)). This is due to the rotation of M caused by the deformation of 
the film, except near the ends where the self-field is operative. The B y near the fixed end of the 



film is nearly equal to the applied field (Fig. 2(d)). This is a consequence of the magnetostatic 
boundary condition, which dictates that the normal component of B will be continuous across a 



boundary. The variation of the magnetic body couple is shown in Fig. 2(e) As the field is applied 
in the transverse direction, in the initial stages the couple is almost uniform throughout the film. 
As the film deforms, the free end gets more aligned with the applied field and the couple intensity 
near the free end is considerably reduced. 

Next, we study the same film but now subjected to a rotating field. The initial field is in the 
x direction and has magnitude -Bomax- Then we linearly increase the field in the y direction and 
decrease it simultaneously in the x direction. By doing so, the magnetic field vector is 'rotated' 
from the x to the y axis. Of course, its magnitude does not remain constant during rotation. 
When the field has rotated, its orientation and magnitude is maintained at a constant value. From 



Figs. 3(a) and 3(b) it can be seen that the deformation is very similar to the previous case (cf. 
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a 




2 


a) 



Figs. 2(a) and 2(b)). At time £ re f (corresponding to instant 5 in Fig. 3(a)), the magnetic forces 
are larger than the elastic forces, hence the film continues to deform for a short time, even when 
the applied magnetic field is kept constant. The contributions to B from various sources is shown 



in Figs. 3(c) and 3(d) Unlike the previous case, the couple intensity always increases at the free 



end (Fig. 3(e)). The effect of the rate of rotation of the magnetic field is shown in Fig. 3(f) We 
decrease the reference time of the field rotation leading to an increased rotation rate. As the rate 
of rotation is increased, the tip reaches its maximum deflection in less time, but its maximum value 
remains the same irrespective of the rate of rotation. 
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(e) 

Figure 2: A horizontal permanently magnetic film subjected to an increasing field in transverse (i.e. vertical) 
direction, (a) The displacement of the free end in time, (b) Deformed geometry: The arrows on the film show the 
direction and magnitude of B and the separate arrows at the right show the applied field Bq. (c) Contributions to 
B x from various sources, (d) Contributions to B y from various sources, (e) Variation of body couple along the film 
at different instances of time. 
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Super-paramagnetic film in a rotating magnetic field. We study the behavior of a SPM film 
subjected to the same rotating applied field as for the permanently magnetic film. The direction 
and magnitude of magnetization in a super-paramagnetic film is not fixed, but changes with the 
applied field and the geometry of the film, given the specific magnetic susceptibility of the film 
(here we have \ x = 4.6 and x y = 0.8). The variation of the tip displacement, the geometry, 
field (both applied and total) and induced couple for a super-paramagnetic film are shown in 



Fig. HI Figs. 4(a) and 4(b) show that for the given susceptibility the deformation of the film is 
similar to that of the permanently magnetic case. As the film deforms, the x component of the 
applied field decreases and the y component increases, and accordingly, so do the components 



of the magnetization and the total field (Figs. 4(c) and 4(d)). It can be observed that the field 



gradients are lower when compared to the permanently magnetic films (see Figs. 



2 c 



2(d), 3(c) 



and 3(d) ). This is due to the fact that the magnetization is not a uniform constant throughout the 
film as in the permanently magnetic case. As the film rotates, unlike the permanently magnetic 
film, the By near the fixed end of the film is not equal to that of the applied field. This is due 
to the fact that field outside the film is enhanced by the normal magnetization of the film. For a 
SPM film, the magnetization is obtained by solving an implicit equation to obtain M (Eqn. 1281) . 
for which it is essential to account for the self-field, generated by the magnetization in the film 
(see Eqn. |27|). Thus, the difference between the applied field B and the magnetic induction B in 



Figs. 4(c) and 4(d) is entirely due to the segment-to-segment interaction. Unlike the permanently 
magnetic film, in the super-paramagnetic film the induced magnetic couple is not uniform, even 
in the initial stages of deformation. As the applied field is rotated, the couple induced at the fixed 
end tends to increase initially, but it gets reduced as the applied field becomes vertical. The couple 
distribution, for a given instant of time, along the film does not monotonically decrease towards 
the free end as in the permanently magnetic film, but attains a maximum along the film. Despite 



this different torque distribution, the deformations are identical (see Figs. 3(b) and 4(b)). 

For the SPM film to deform, the material susceptibility tensor x needs to be anisotropic. The 
induced magnetic couple depends on the difference between \ x and \ y and on the magnetic field, 



14 




Figure 3: A horizontal permanently magnetic film subjected to a rotating field, (a) The displacement of the free 
end in time, (b) Deformed geometry: The arrows on the film show the direction and magnitude of B and the 
separate arrows show the applied field, (c) Contributions to B x from various sources, (d) Contributions to B y from 
various sources, (e) Variation of body couple along the film at different instances of time, (f) Displacement of the 
free end vs time for various loading rates t re f. 
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(e) 



Figure 4: A horizontal super-paramagnetic film subjected to a rotating field in transverse direction, (a) The 
displacement of the free end in time, (b) Deformed geometry: The arrows on the film show the direction and 
magnitude of B and the separate arrows show the applied field, (c) Contributions to B x from various sources, (d) 
Contributions to B y from various sources, (e) Variation of body couple along the film at different instances of time. 
To normalize the couple distribution in Fig. |4(e) the magnetization is used of the permanently magnetic film (see 
Fig. [3(e) I. 
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which itself depends on the magnitude of the susceptibility. The effect of the change in Xx, with 



constant Xyi on the deflection of the free end is shown in Fig. 5(a) As Xx decreases, the induced 
magnetization decreases, and so does the magnetic couple, resulting in less deflection of the film. 
Also, for small values of Xx (< 2.0), even when the applied magnetic field is kept constant after 
t ie {, the magnetic forces due to the induced couples are smaller than the elastic forces. Hence, the 
film returns back to its initial position. The deflection of the free end with varying Xy an d Xx 



but keeping their difference constant is shown in Fig. 5(b) As the magnitude of the susceptibility 
increases, the deflection of the film decreases. Hence, for a SPM film both the absolute value of 
the susceptibility tensor as well as its anisotropy determine the deformation of the film. 



The effect of the rate of rotation of the magnetic field is shown in Fig. 5(c) We reduce the 
time in which the applied field is rotated from horizontal to vertical (t re f) from 1 ms to 0.4 ms. 
As a result the rotation rate of the applied magnetic field is increased. This study is done for two 
different cases of anisotropy: Xx = 4.6, Xy — 0.8 and Xx = 2.0, Xy — 0.8. When Xx = 4.6 and 
X y = 0.8, in the steady state, the deflection of the free end attains a value independent of the rate 
of rotation of the applied magnetic field. The only effect of the increase of the rate of rotation, 
the free end attains its maximum displacement in less time. When Xx — 2.0 and Xy = 0.8, 
however, the maximum deflection of the free end decreases when the rotation rate is increased. 
This behaviour is typical for the SPM film. The magnetic torque distribution depends sensitively 
on the orientation of the film with respect to the field direction. If the field is rotated too fast, 
the film lags behind due to the opposing fluid drag. A field oriented perpendicular to an almost 
undeformed film, will not induce a pronounced torque and the film returns back to its undeformed 
configuration by elastic recovery. This is the case when the susceptibility is unfavourable as for 



(Xx, Xy) = (2-0, 0.8) in Figs. 5(a) and 5(c) and for (x x , Xy) = (5-0, 2.8) in Fig. 5(b) 



3.2. Dimensional analysis 

To identify the dimensionless parameters that govern the behavior of the system, we start from 
the virtual work equation (see Section [27TT) . neglecting the axial deformations: 



,d 2 vd 2 5v , f .d 2 v . , f Ar d5v . , f _ dv 

d^^ dx + / p A w 5vdx = / N ^ Adx - / c «di l 
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Figure 5: A horizontal super-paramagnetic film subjected to a rotating field: (a) Displacement of the free end vs 
time for various \ x {\ y = 0.8). (b) Displacement of the free end vs time for various \ x and \ y , with Xx — Xy = 2-2. 
(c) Displacement of the free end vs time for various t Te {. 
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where the first term represents the virtual elastic work done by the internal moments, the second 
term represents the virtual work done by the inertial forces, the third term represents the virtual 
work done by the magnetic couple and the last term represents the work done by the fluid drag 
forces. This is valid for any segment dx, and hence, 

rnT d 2 vd 2 Sv A d 2 v . , r 85v . dv . , 

We introduce the dimensionless variables V, T and X, such that v = VL, x = XL and t = Tt, 
where L is a characteristic length (taken to be the length of the film) and r is a characteristic 
time. Substitution yields 

Ebh 3 d 2 Vd 2 5V pbhL 2 ^d 2 V C v L 2 bdV „ r ^ 17 d5V 

^imlm + — 5v lm = -^rw sv + N ' hb l)x- (29) 

from which the elastic (ff^r), the inertial ( pb ^f ), the viscous ( Cy ^ h b ) and the magnetic (N z hb) 
terms can be easily identified. By normalising with the elastic term, we get 

1 fd 2 Vd 2 5V\ T fd 2 V^\ „ (d5V\ „ (dV„\ 

With the three governing dimensionless numbers being defined as, I n = £h 2 the ratio of 
inertial to elastic force, M n = the ratio of magnetic to elastic force and F n = the ratio 
of fluid to elastic force. 

Before proceeding, we identify the origin of the magnetic couple N z for the two magnetic 
material systems under consideration. For permanently magnetic materials, having a remanent 
magnetization along the axial direction of the film:(M x , M y ) = (M x , 0), 

N z = M x B y - M y B x = M x B y = M x B f(9), (31) 

where 8 is the film orientation and B is the amplitude of the applied magnetic field. For a SPM 
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film, 



N z =M x B y - M y B x = n H x H y (xx - Xy) 

Xl Xx )9x (0)H o f y (x x , Xy)9y(0)(Xx ~ Xy) (32) 

B 2 - 

=—h(Xx,Xx,0), 

in which h groups the dimensionless dependence of the torque on the susceptibilities Xx-, Xyi 

the 

film orientation 9. It is to be noted that the body couple for the SPM film is proportional to the 
square of the applied field while in the permanently magnetic case, it is linearly proportional. As 
a result, the magneto-elastic number becomes M n = m ^° 2 l2 for the permanently magnetic case 
and M n = ^ o °^ h2 for the SPM case. 

Fig. El shows the normalised tip displacement as a function of time for permanently magnetic 
film subject to a rotating field, as studied in Section 13.11 

The transverse displacement is plotted with varying each of the non-dimensional number keep- 
ing others constant. The reference cases (M n = M nre f, I n = I nre { and F n = F nre f) are the same as 
those shown in the fundamental loading situations. 

It is to be noted from Figs. El El and [TJ that the steady state of the film deformation is principally 
governed by the ratio of magnetic to elastic forces. With increase of the magneto-elastic number 



the film deformation increases (Figs. 6(a), 7(a)). In the case of SPM, for low magneto-elastic 



number the elastic forces are much larger than the magnetic forces. Hence the film fails to deform 
after a certain extent and returns to its initial position. 

Changing the fluid and inertia numbers changes the dynamic properties of the system ( damping 
and time period). Increasing the fluid number changes the system from under-damped to an over- 
damped state. Increasing the inertia number changes the time period of oscillation by increasing 
its period and increases the amplitude of oscillation. When I n is small (/ = 0.1/ re f), the film shows 
a near quasi-static deformation behavior. 

The increase of the inertia number makes the film heavier, hence the film moves slowly, showing 
an increase in its time period. As the film has higher inertia number, the inertia forces will be high, 
due to which the film undergoes large amplitudes oscillation about the steady state. Another way 
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Figure 6: Transverse tip displacement of a permanently magnetic film with variation of non-dimensional numbers, 
(a) M n , (b) F n , (c) I n , 

to look at it is, increase of the inertia number decreases the elastic forces (I n =Inertia force/elastic 
force) resisting the deformation of the film, thus making the film floppy, hence its amplitude of 
oscillation about the steady state is large. The increase of fluid number increases the drag force, 
hence the deformation of the film in the transient stage is less, thus taking a long time for the film 
to reach the steady state. These two facts (larger inertia number leads to large deformation and 
large fluid number leads to less deformation) will be of help in understanding the behavior of the 
asymmetric configurations. 

As can be seen in Figs. El and [HI three modes of deformation can be identified. The 
under-damped motion of the film, where the film oscillates around the steady state before reaching 
it, the over-damped motion, where the film takes a long time to reach the steady state and the 



quasi-static motion of the film where the dynamic effects are not seen. It is seen from Figs. 6(a) 
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and 7(a) the amplitude of film deformation in the steady state is governed by the magneto-elastic 
number (M n ). The inertia number and the fluid number determine the way in which the steady 
state is reached (under-damped, over-damped or quasi-static). 

The motion of the film can be quantified with the steady state deflection, the time to reach the 
steady state and the amplitude of oscillation. The latter two are independent of the steady state 
deflection and can be used to quantify the modes of deformation according to Table CD 

We now extend the range of fluid and inertia number so as to include all possible materials, 
fluids and geometries. This is done for the case of permanently magnetic materials. The actuator 
properties chosen cover the length of the film from 100 nm to 1 mm, film densities from 1000 to 
8000 kg/m 3 , elastic moduli from 1 MPa to 200 GPa and fluid viscosities from one tenth to ten 
times that of water. As a result the fluid and and inertia number are scanned by 8 orders of 
magnitude, for a given magneto-elastic number. For a given magneto-elastic number the time to 
reach the steady state (t ss ) and the amplitude of oscillation are plotted against varying inertia and 
fluid number in Fig. [9j 

From the contour plots of the amplitude of oscillation, the region in which the film shows 
oscillatory behavior before it reaches a steady state can be identified. In this region, the plots of 
the time to reach steady state show that the film reaches the steady state after the field has stopped 
rotating. This under-darned behavior occurs for large inertia number and low fluid number. 

From the contour plots of the time to reach the steady state, a region can be identified when 
the film reaches the steady state as quickly as possible with no oscillation. Absence of oscillation 
implies that the inertia effects are negligible and because the film reaches the steady state as 
quickly as possible, the time lag due to viscous effects are absent and hence the viscous forces 
are negligible. The film deforms because of the competition between elastic and magnetic forces 
and the dynamic effects are negligible. Such a deformation behavior is named as "quasi-static". 
Quasi-static behavior is observed when the inertia and fluid numbers are small. 

From the contour plots of amplitude of oscillation and the time to reach the steady state, 
another region can be identified in which the film shows no oscillatory behavior, but reaches the 
steady state after the field has stopped rotating. As the film exhibits no oscillations, the inertia 
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Figure 8: Variation of transverse displacement of a permanently magnetic film with magneto-elastic number for 
different systems identified by the symbols shown in Fig. 9(a) 
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Figure 9: A/„=0.4688. (a) Contours of amplitude, |(b)| Contours of amplitude and time to reach steady state. 



effects are negligible in this region. The fluid number in this region leads to large viscous forces, 
which make the film to deform in an over-damped manner. 

The transitions between the zones happens as follows. If the system is initially in the over- 
damped region, reducing the fluid number leads to the following consequences: At large inertia 
numbers, the system will start to show oscillations about the steady state and enter the under- 
damped region. At low inertia numbers, the t ss gets reduced and at a certain fluid number the 
film reaches steady state at t re f. Further reduction of the fluid number does not change the t ss and 
brings the system in to the quasi-static region. 

Decrease of the inertia number also leads to two kind of transitions: At large fluid numbers, it 
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Table 1: Deformation modes quantified 



Deformation mode 


time to reach 


steady state t ss 


amplitude 


under-damped 


> 


tref 


> 


over-damped 


> 


tref 





quasi-static 











Fluid number 

Figure 10: Schematic representation of the deformation behavior of the film at different inertia and fluid numbers. 
Systems along the solid lines have constant shape and the systems along the dashed lines have constant size. 

reduces the amplitude of oscillation, and, as the viscous effects are large for large fluid numbers, 
enables the system to enter the over-damped region, where no oscillations are observed. At low 
fluid numbers, the decrease of inertia number decreases the amplitude of oscillation, and, as the 
viscous effects are small, the system is shifted to the quasi-static region, where the dynamic effects 
are negligible. 

The deformation modes can be nicely summarized in terms of the dimensionless numbers as 
shown in the schematic picture shown in Fig. [101 The transition from over-damped to quasi-static 
region is obtained from the plots of steady state time (Fig. [9]). The transition from the under- 
damped to the other two regions is obtained from the plots of amplitude of oscillation. For a given 
system, this diagram can quantify its the behavior at any magneto-elasic number. 

The effect of size is investigated now. Keeping the ratio L/h fixed, their magnitudes are varied 
and the corresponding displacements are plotted for permanently magnetic and SPM case (Fig. fTTj) . 
With reference to the non-dimensional numbers, keeping the ratio L/h fixed and varying their size 
does not change the magneto-elastic number, but has noticeable effect on the fluid number and 
the inertial number. When the size of the film is small, it shows a near quasi-static behavior, 
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Figure 11: Transverse tip displacement of film, with same L/h ratio while their magnitudes are varied, (a) Perma- 
nently magnetic film, (b) Super-paramagnetic film. 

L = 10/zm in Fig. [UJ Which enables us to infer that, as the size of the film decreases the inertia of 
the film can be neglected. This is similar to the fluid dynamic counterpart, where at small length 
scales the inertia and viscous effects are negligible. 

Referring to the schematic diagram in Fig. [101 the solid lines on the plot shows the effect of size 
(for a given modulus, density, fluid, characteristic time) on the actuator system. The solid lines 
are the lines of constant shape and as we decrease the size (move in the direction of the arrow), 
the fluid and inertia number are reduced. Thus the system moves towards the quasi-static region 
when the size of the actuator is reduced. It can be shown that these lines of constant shape have 
a slope of 2. 

The dashed lines on the schematic diagram represent the lines of constant length when the 
aspect ratio of the film is varied (for a given modulus, density, fluid, characteristic time). As the 
aspect ratio of the film is decreased (move in the direction of the arrow) it becomes more stiff and 
the elastic forces in the film increase, thus reducing the inertia and fluid number. Thus, taking 
the film to the quasi-static region. It can be shown that these lines of constant size have a slope 



3.3. Asymmetric motion 

As mentioned before, the fluid propulsion in micro-channels takes place at low Reynolds num- 
ber. The actuating member, therefore, should move in an asymmetric manner to effectively propel 
fluid. In this Section we discuss several configurations that are able to do so. In all the configura- 



of 2/3. 
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tions the applied magnetic field is uniform in space, but its magnitude and direction are varied in 
time. By tuning the applied field, the initial geometry of the film and its magnetic nature (perma- 
nently magnetic or super-paramagnetic), we have identified four configurations that mimic ciliary 
motion. For the results presented, the thickness of the film is taken to be 2 fim and the elastic 
modulus to be 1 MPa. The drag coefficients are calibrated against computational fluid dynamic 
(CFD) simulations, details of which will be discussed in Section l3~4l 



1. Partly magnetic film with cracks. The natural cilium is found to have a varying stiffness 
in the effective and recovery stroke [5|. To use this concept we need to have the film to possess a 
large bending stiffness in the effective stroke while pushing the fluid and to possess a low stiffness 
during the recovery stroke. This can be achieved by introducing cracks in one side of the film, 
while only a part of the film is magnetic. The film is straight initially, is attached at the left end 
and has cracks of size 0.75 /im at the bottom. By magnetizing onlya part of the film, the film is 



expected to behave like a flexible oar (as also mentioned by Purcell |12|). Only 20% of the film, the 
end near to the fixed part, is magnetic. The assumed remnant magnetization is 15 kA/m, with the 
magnetization vector pointing from the fixed end to the free end. The drag coefficients used are 
C x = 30 Ns/m 3 and C y = 60 Ns/m 3 . The applied magnetic field is increased linearly to 145 mT in 
the y direction in 0.6 ms, then rotated by 90° in the next 1.2 ms and finally reduced to zero in the 
next 0.2 ms. These system properties lead to a magneto-elastic number of 5.43 or 1.08 (0.2 x 5.43), 
a fluid number of 0.416 and an inertia number of 0.012. The movement of the film under the action 
of the applied magnetic field is shown in Fig. [121 When the external magnetic field is applied, the 
magnetic couples act on the magnetized portion of the film in a counter-clockwise manner, thus 
rotating the film about the fixed end. Now the drag forces are acting on the top part of the film, 
which close the cracks, making the film stiff. 



When the applied field is switched off (Fig. 12(e) ), the film will recover elastically and the drag 
forces act on the bottom part of the film which open the cracks making the film compliant. Such 
an interaction of magnetic couples, elastic forces and drag forces results in an asymmetric motion, 
as is clear from Fig. [121 



27 



(a) t = 0.2ms (b) t = 0.4ms (c) t = 1.0ms (d) t = 1.9ms 




(e) t = 2.1ms (f) t = 3.0ms (g) t = 3.5ms (h) t = 6.5ms 




(i) Trajectory of the free end 

Figure 12: Film with cracks with only a part (20%) magnetized. The dashed line shows the initial position of the 
film. The arrow shows the direction of the applied field. 



2. Buckling of a straight magnetic film. A straight horizontal magnetic film with a slight 
perturbation is used to get the desired asymmetric motion. The film is assumed to have a uniform 
magnetization with the magnetization vector pointing along the film length, from the fixed end at 
the left to the free end at the right. The remnant magnetization of the film is taken to be 15 kA/m. 
The length of the film is 100 /iin. The drag coefficients used are C x = 30 Ns/m 3 and C y = 60 
Ns/m 3 . The external field is applied as follows: A field of 30 mT is applied in the negative x 
direction from t — to t — 1 ms and the field is reduced to zero in the next 0.2 ms. These system 
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properties lead to a magneto-elastic number of 1.125, a fluid number of 0.75 and an inertia number 
of 0.04. Initially, the magnetization and the applied field are parallel, but with opposite sign, so 
that the magnetic couple is zero. However, with any perturbation of the film, the equilibrium state 
becomes unstable and the film will buckle away from the straight configuration. By assuming a 
uniform magnetization in the film and neglecting drag forces, the critical field can be calculated 
(see Appendix [B]) . When it buckles the film tries to curl such that the magnetization in the film 





(i) Trajectory of the free end 

Figure 13: Movement of a perturbed film with the applied magnetic field in opposite direction to the magnetization 
of the film. The dashed line shows the initial position of the film. The arrow shows the direction of the applied 
field. 



is aligned with the applied field (Fig. 13(a) 13(d)). When the applied field is removed, the film 



returns back to its initial position by elastic recovery. For this configuration, the effective fluid 
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propulsion will take place when the film recovers elastically. Whether the film will buckle up or 
down depends critically on the sign of the initial imperfection (see Appendix [B]) . 



3. A curled magnetic film. A curled film with remnant magnetization is subjected to a uniform 



magnetic field. The initial geometry of the film is shown in Fig. 14(a) The left edge is the clamped 
end. The direction of the magnetization is along the film with the magnetization vector pointing 
from the clamped end to the free end. The remnant magnetization of the film is taken to be 15 
kA/m. An external field of magnitude 9 mT is applied at 225° to the x axis from t = ms to 
t — 1 ms and then linearly reduced to zero in the next 0.2 ms. The drag coefficients used are 
C x = 5 Ns/m 3 and C y = 60 Ns/m 3 . The radius of curvature of the film is 100 fim. These system 
properties lead to a magneto-elastic number of 0.831, a fluid number of 4.55 and an inertia number 
of 0.24. The physical mechanism of asymmetry is akin to the previous case, except that the film 
is made to buckle in a predetermined way by a curled geometry, which turns out to exhibit a 
large asymmetry in motion. The propulsive action in the effective stroke takes place during elastic 
recovery. 

In the portion near the fixed end, the magnetic couple acts in a clockwise sense and in the 
portion near the free end, in a counter-clockwise sense. As a result, these couples tend to bend the 
film, such that the curvature of the film increases, bringing the ends of the film closer together (this 
behavior is similar to that of a natural cilium). Here, the asymmetry is larger than the previous 



case as can be observed from Figs. 13 (i) and 14(i) 



4. Super-paramagnetic film. A super-paramagnetic film which is anisotropic in its magnetic 
susceptibility and which is tapered along its length is subjected to a rotating magnetic field. The 
assumed susceptibilities are 4.6 and 0.8 in the tangential and normal directions, respectively. The 
thickness of the film varies linearly along its length, being 2 /xm at the left (attached) end and 
decreasing to 1/xm at the right end. A magnetic field of 31.5 mT is rotated from 0° to 180° in 
t = 10 ms and then kept constant for the rest of the time. The drag coefficients used are C x = 30 
Ns/m 3 and C y = 60 Ns/m 3 . These system properties lead to a magneto-elastic number of 1.96, 
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(a) t = 0.0ms 



(b) t = 0.1ms 





(e) t = 1.1ms 



(f) t = 1.5ms 



(c) t = 0.2ms 




(g) i = 2.0ms 




(d) t = 0.3ms 



(h) t = 6.0ms 




(i) Trajectory of the free end 
Figure 14: Movement of a curled permanently magnetic film. The arrow shows the direction of applied field. 

a fluid number of 0.075 and an inertia number of 4.0 x 10 -4 . When the rotating field is applied, 
the free end portion of the beam is rotated through 135° and a U bend is formed near the fixed 



end (Fig. 15(e) and 15(f)). It is to be noted that the portion of the beam near the free end is 



nearly straight. This is because the magnetization of the film in this region is almost aligned to 
the applied field (which is evident from the magnetization vectors shown on the film in Fig. [T5|) . 
As a result, the normal component of the field is low, and hence the moment is less as well. In 
this situation it is in the part of the film near the U bend, where the magnetic couple distribution 
(which tends the film to bend) balances the elastic forces (which tends the film to become straight 
again), hence " freezing- in" the bent shape. When compared with other parts of the beam, at the 
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U bend portion of film the magnetic couples are large. On the portion of the film between the free 
end and the U bend, anticlockwise moments are acting. On the portion between the fixed end and 



the U bend, clockwise moments are acting (Figs. 15(e) , 15(f) , 15(g) ). Due to the tapered nature of 



the film, the clockwise moments are larger than the anticlockwise moments. Under the influence 
of such a system of moments, the film becomes more curved, the end to end distance decreases and 
the film recovers. As the beam recovers elastically the U bend propagates to the free end of the 
beam. It is to be noted that the film recovers in the presence of magnetic forces, i.e. the recovery 
is not an elastic one, but controlled by magnetic forces, keeping the film low. This phenomena can 
be exploited to provide a large asymmetry in motion of the film during the forward and return 
stroke. A larger field will result in larger curvature of the U bend, forcing the film to stay lower, 
enhancing the efficiency of the return stroke. 

3.4- Fluid propulsion 

In the magneto-mechanical model, the presence of the fluid is incorporated through drag rela- 
tions that account for the surface tractions as a function of the velocity. The specific values for 
the drag coefficients used in the previous Section have been calibrated through a comparison with 
a fully-coupled solid-fluid model. This model and the calibration procedure will be discussed in 
this Section. Moreover, we will also investigate how efficient the four asymmetric motions are in 
actually propelling fluid. We will demonstrate that the swept area by the film tip is one-to-one 
related to the fluid propelled. 

The Lagrangian solid dynamics model used to study the magneto-mechanical behavior of the 
films (Section [2]) is coupled to an Eulerian fluid dynamics code based on the method detailed in 
Fluid inertia is neglected, so the fluid model effectively solves the Stokes equations. The 
explicit coupling between the two domains is established through Lagrange multipliers. Input to 
the fluid dynamics model are the positions and velocities of the film at all times which result in a 
full velocity field in the fluid. We calculate the drag forces on the film as tractions via the stress 
tensor in the fluid. The traction distribution is subsequently imposed as surface tractions in the 
solid dynamics model through Eqn. [13j 

To calibrate the drag coefficients, coupled solid-fluid simulations are performed using a periodic 
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(a) t = 1.0ms (b) t = 3.0ms (c) t = 5.0ms (d) t = 6.0ms 




(e) t = 7.0ms (f) t = 8.0ms (g) t = 8.5ms (h) t = 10.0ms 




(i) Trajectory of the free end 

Figure 15: Super-paramagnetic film in a rotating magnetic field. The dashed line shows the initial position of the 
film. The big arrow shows the direction of the applied magnetic field. The small arrows show the magnetization 
along the film. 

arrangement of cilia in a micro- fluidic channel. The dimensions of the unit-cell analyzed are 400/im 
in width (horizontal) and 500/zm in height. No-slip boundary conditions are applied at the top 
and bottom boundaries of the channel and periodic boundary conditions at the left and right ends 
of the unit-cell. The viscosity of the fluid is taken to be that of water (1 x 10~ 3 Pas). The film is 
placed at the center of the channel to avoid any interaction with the boundary. From the coupled 
solid-fluid simulations we get the trajectory of the free end. We then perform simulations using the 
magneto- mechanical model (Section [2} with the assumption that the fluid exerts tractions that are 
proportional to the velocity, on the film as mentioned in Section 12.11 We vary the drag coefficients 
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to match the trajectory of the free end obtained from the coupled solid-fluid simulations. The 
results are depicted in Fig. [161 clearly showing that the agreement between the uncoupled and 
fluid-coupled trajectories is very good. The drag coefficients obtained here have been used for the 
simulations presented in Section l3~3l 
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(a) Partly magnetic film with cracks. 




(b) Buckling of a straight magnetic film. 
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(c) Curled magnetic film. 




(d) Super-paramagnetic film 



Figure 16: Trajectory of the free end for various cases obtained when calibrating the drag coeffients used in the 
solid-drag simulations of Section 13.31 with the fully coupled solid-fluid simulations. For the magnetic film with 
cracks, the crack length is taken to be 0.1/xm; the parameters of the three other configurations coincide with Section 
3^1 



At low Reynolds numbers, asymmetric motion is required to propel the fluid. To study how 
much fluid is propelled for a given asymmetric motion, we studied two configurations: the curled 
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permanently magnetic film (Fig. [H]) and the SPM film (Fig. [15!) . As a measure for the asymmetry, 
we compute the area swept by the free end of the film. To vary this area we vary the magnitude of 
the magnetic field applied and perform coupled solid-fluid simulations. The corresponding swept 
area and the total flow across the channel in one cycle is computed. The domain for which the 
calculations are performed, is the same as that used for the calibration of the drag coefficients, 
discussed above. The dependence of fluid propelled (the area flow per cycle) on the swept area 
is shown, for both the cases, in Fig. [I7J The cycle times for the two cases are 15 and 10 ms for 
the curled magnetic and the super-paramagnetic films, respectively. The volumetric flowrate can 
be calculated from Fig. [17] by multiplying the area flow per cycle with the out-of-plane dimension 
of the channel divided by the cycle time. The flow rate is seen to vary linearly with the swept 
area. This suggests that the swept area can be used as a measure of effectiveness of the cilium, 
representing the fluid volume displaced. In the following Section, we will use this to investigate 
the efficiency of fluid propulsion as a function of the system parameters. 




3 3.3 3.6 3.9 4.2 2 3 4 5 6 

Swept area (1 CT 3 mm 2 ) Swept area (1 CT 3 mm 2 ) 

(a) Curled magnetic film. (b) Super-paramagnetic film. 

Figure 17: Relation between swept area and fluid propelled. 



4. Efficiency survey and discussion 

A parametric study of the curled magnetic film and the super-paramagnetic film is performed 
next. The parameters that govern the response of the actuator system are the inertia number, 
the fluid number and the magneto-elastic number. The performance of the system is quantified in 
terms of the swept area (as a measure of the asymmetry and associated fluid flow) and the time 
taken by the film to reach its initial position. 
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The actuator system based on its material properties, its geometry and the fluid to be propelled 
can have different fluid and inertia numbers. Fixing them establishes the system. Different set of 
these numbers refer to different systems. For a set of these systems the magneto-elastic number is 
varied and the swept area and the time to reach the initial position is noted. 

The ingredients in the dimensionless parameters are as follows: L is taken to be the length of 
the film, h is the thickness of the film, the length is varied from 100 to 1000 /im and the aspect 
ratio {L/h) is varied from 10 to 200. The density p of the film is varied from 600 to 8000 kg/m 3 , 
the drag coefficient C y is varied from 6 to 600 Ns/m 3 , while the ratio C y /C x is kept the same, and 
the time scale r is taken to be the time in which the load is applied. For a magnetic film it is the 
time during which the magnetic field is applied and for the SPM it is the time during which the 
magnetic field is rotated. 

The area swept by the free end of the film for different systems (varying inertia number and 



fluid number) keeping the magneto-elastic number fixed is shown in Figs. 18(a) and 19(a), for the 
magnetic film and SPM film respectively. The area is normalized with the maximum area that can 
be swept by a given system (maximum area that can be swept is 7rL 2 /2). This plot shows that 
for a given magneto-elastic number, the swept area decreases with an increase of fluid number. 
As the fluid number increases, the drag force opposing the motion of the film increases, hence the 
deformation of the film is less, thus sweeping lesser area. 



The magnetic number needed to sweep a given area for different systems is shown in Figs. 18(b) 



and 19(b) It can be seen that the same area can be swept by all the systems. As the fluid number 
increases, we need a large magneto-elastic number to sweep a given area , for a given inertia 
number. When the fluid number is kept constant, we need a larger magneto-elastic number to 
sweep the same area at low inertia numbers. For the curled magnetic film, dependency of the 
swept area of the film is large at small fluid numbers and this dependency decreases at high fluid 
numbers. This is due to the fact that the viscous forces dominate over the inertial forces at high 
fluid numbers. 



When sweeping an area of Fig. 18(b) and 19(b), time taken by different systems to return 



to their initial position is shown in Figs. 18(c) and 19(c), respectively. It can be seen that it 
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takes longer for the film to return to its initial position when the fluid number is increased and is 
independent of the inertia number. But in the case of SPM film, because of the whip like return 
stroke, some dynamic effects are observable at larger inertia numbers, hence these systems take 
longer time to return to their initial position. 

The lines of constant magneto-elastic number for different systems which sweep the same area 



are shown schematically in Fig. 18(d) and 19(d) for the magnetic film and the super-paramagnetic 



film, respectively. These figures nicely summarize the dependency of swept area on the non- 
dimensional parameters. The arrow on broken line points in the direction of increasing magneto- 
elastic number, which means that to get the same swept area we need higher magneto-elastic 
numbers at high fluid numbers. The lines of constant magneto-elastic number have a slope, this 
shows the dependency on the inertia number, i. e. as the inertia number increases, to sweep the 



same area we need relatively less magneto-elastic number. For the curled magnetic film (Fig. 18(d) ), 
the change in slope of the constant magneto-elastic number implies the decrease in dependency on 
the inertia number at high fluid numbers. For the curled magnetic film, the upper bound for the 



Inertial number, shown in Fig. 18(d) has a slope of 2.0 The lines of constant shape (L/h constant) 
and constant size (L constant) are also shown in the schematic pictures. It can be seen that 
reducing either the size or the aspect ratio decreases the fluid number. Hence to sweep a given 
area the film needs lower magnetic number when either the aspect ratio of the film or its size is 
decreased. 

5. Conclusion 

This work has addressed the numerical design of magnetically driven actuators for fluid propul- 
sion in micro-fluidic channels. The key characteristics of such actuators, is that they have to move 
in an asymmetric manner, due to the low Reynolds number. Such actuators are proposed to be 
made of polymer films with embedded magnetic particles which can be manipulated by an exter- 
nally applied field. Four possible configurations have been identified which result in an asymmetric 
motion of the actuators. In the configuration based on permanently magnetic films the effective 
stroke is slower compared to the recovery stroke. This is due to the fact that the effective stroke 
is driven by the elastic forces, while the recovery is due to the applied magnetic field. This nature 
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Figure 18: Parametric study of a curled permanently magnetic film, (a) Contours of swept area for M n — 30.0, |(b)| 
Contours of M n needed to sweep an area of 0.2, (c) Time taken by the film to return to its initial position when the 



tip sweeps an area of 0.2, and |(d)| Schematic representation of the lines of constant magneto-elastic number. The 
arrow on the broken, dashed and solid lines represent the direction of increasing magneto-elastic number, decreasing 
aspect ratio and decreasing size, respectively. 
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Figure 19: Parametric study of a curled permanently magnetic film, (a) Contours of swept area for M n — 30.0, |(b)| 
Contours of M n needed to sweep an area of 0.45, (c) Time taken by the film to return to its initial position when the 
tip sweeps an area of 0.45 and | (d) | Schematic representation of the lines of constant magneto-elastic number. The 
arrow on the broken, dashed and solid lines represent the direction of increasing magneto-elastic number, decreasing 
aspect ratio and decreasing size, respectively. 
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of slow effective stroke and fast recovery stroke has also been reported by Kim and Netz 9j, in 
their study of the propulsion efficiency of a periodically beating elastic filaments anchored to a 
substrate. An important result from our work is that the flow across the channel varies linearly 
with the area swept by the free end of the film, promoting the swept area as an efficiency parameter 
of the actuator system. The permanently magnetic film is very sensitive to the applied field at 
large values of magnetization. The SPM film is equally sensitive to the applied field at all values 
of susceptibility. 
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Appendices 

A. Interpolation Functions 



N u = [N x N 2 0] 
N v = [o E x H 2 H 3 H A ] , 
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where 



Ni =1 - 



x 

T 
1 



Hx =^(2x 3 -3x 2 l + l 3 ), 



Ho 



/n/ 3 



yi - 2x 2 t 2 + xi 6 ) 



H 3 =-{-2x 3 + 3x 2 l), 



Ha 



LI 3 



2r2\ 



and I is the length of the element and Iq. is a reference length. 

B. Magnetic Buckling Analysis 

In this Section we study the buckling behavior of a straight permanently magnetic film. When 
the applied field is opposite to the magnetization and the film is straight, no couple is induced for 
the motion of the film (see Fig. |20{) . But, if the film is slightly perturbed, couples will be acting 
on the film. These couples will increase when the film is deflected further from the initial, straight 
configuration. Clearly, the straight configuration is an unstable equilibrium state, so that buckling 
will occur above a critical value of the magnetic field. This situation is similar to a cantilever 
with an end compressive load. The cantilever buckles because of the moment of compressive load 
and the film buckles because of the magnetic couple. The total moment of the body couple on an 

M 




Figure 20: Buckling of a film under a magnetic field 
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infinitesimal element of length dx at x is 



dN = MAB sin 6dx, 



where 9 is the angle between the beam axis and the horizontal. The total moment at x 



x rx fx [X j 

N(x) = I dN x dx= / MAB sin 9dx = / MABOdx = / MAB—dx 
'o Jo Jo Jo 

fX 



dx 



[ MABdv = MAB (v(x) - v ) 
Jo 



The equilibrium equation of the beam is 

.d 2 v (x 

dx 2 
d 2 v(x) 

dx 2 

d 2 v(x) 



El-^r- = -MAB(v(x)-v ) 



-k (v(x) - D ) 



dx 2 



+ k v(x) = k vq, 



where, k 2 = and EI is the bending stiffness of the beam. The solution is 



v(x) = A cos kx + B sin kx + vq. 



Using the boundary condition at the clamped end 



we get, 

v(x) = v (l — cos k(l — x)). 
The displacement at x = is v , which yields 

v o — Vq(1 — cos kl), 
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so that the buckling criterion can be written as 

EIn 2 n 2 , . 

MB = _-. (30) 

For a given M, the critical B can be found from Eqn. [B] Note that this kind of behavior cannot 
be observed in a super-paramagnetic film. 
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